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Abstract. Chemical signals in animals stimulate different links and relationships. In lizards, chemical signals play an essen- 
tial role in reproductive and social communication. Although lizard species inhabit tropical ecosystems in great diversity, 
the chemical composition of the femoral gland secretions that tropical lizards use for communication is poorly known. In 
this work, we used gas chromatography-mass spectrometry (GC-MS) to identify and compare the main compounds con- 
tained in the femoral gland secretions of the males of three species of arboreal lizards, Hemidactylus frenatus, Sceloporus 
melanorhinus, and Urosaurus bicarinatus; and of two species of terrestrial lizards, Aspidoscelis lineattissimus and Sceloporus 
pyrocephalus. We identified twelve lipophilic compounds contained, corresponding to fatty acid esters (~58%), carboxylic 
fatty acids (~25%), and other compounds (~17%). Carboxylic acids, mainly with esters content, revealed significant differ- 
ences between lizard species. The potential role of these chemical compounds to social and sexual behaviours is discussed. 
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Introduction 


Chemical signals are pivotal information resources for 
communication in several animal groups (MULLER- 
SCHWARZE 2006, WYATT 2014). The use of chemical signals 
in both inter- and intraspecific interactions are of major in- 
terest to the study of animal behaviour (MARTÍN & LOPEZ 
2014, WYATT 2017). In reptiles, for example, some of these 
semiochemicals contribute to the control of social and re- 
productive behaviours (HoucK 2009, MASON & PARKER 
2010, MARTIN & LOPEZ 2011). 

In recent years, lizards have become a popular mod- 
el group to study the composition and function of verte- 
brate chemical signals (e.g., BAECKENS 2019). Lizards pro- 
duce and release several chemical compounds through the 
skin or specialized body glands (WELDON et al. 2008). For 
example, many species have femoral or precloacal glands 
(MARTÍN & LOPEZ 2011, KHANNOON 2012, VALDECANTOS 
et al. 2014). These glands produce holocrine secretions that 
are released through epidermal pores, and are constitut- 
ed mainly of lipids and proteins (MARTÍN & LÓPEZ 2015, 
MANGIACOTTI et al. 2019). Both male and female lizards 
have femoral (or the analogous precloacal) pores, but in 
many species, secretions are only observed in males, which 


have more abundant secretions during the mating season 
(MARTIN & LOPEZ 2015). 

Although tropical ecosystems are home to a great di- 
versity of lizard species (PIANKA & VITT 2003), the semio- 
chemicals that many lizards use for communication are 
understudied compared to lacertid species from Mediter- 
ranean environments (GARCÍA-ROA et al. 2016, BAECKENS 
et al. 2018). 

In the present study, we identify the lipophilic fraction 
of chemical compounds in the femoral secretions of males 
of three arboreal and two terrestrial lizard species from a 
tropical area of western Mexico, and evaluate the existence 
of interspecific differences in the relative abundance of 
each group of chemical compounds. We discuss our results 
with previously reported chemical compounds in other 
species from similar or different environments. 


Materials and methods 
Study area and lizard species 


During fieldwork carried out between May and June 
2017, coinciding with lizards’ mating seasons (RAM{REZ- 
BAUTISTA et al. 1995, 2000, 2006, RAMfREZ-BAUTISTA & 
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OLVERA-BECERRIL 2004, Diaz et al. 2017), we captured by 
loop adult males of five species in the vicinity of Maruata 
Bay in the Michoacan State on the west coast of Mexico 
(18.271944° N, -103.355556° W). This area is characterized 
by a dry tropical climate with an average annual tempera- 
ture of 28°C and a long annual dry season from Novem- 
ber through June (ALVARADO et al. 1995). The dominant 
vegetation is dry deciduous forest (RZEDOWSKI 1978, AL- 
VARADO et al. 1995). There are also patches of xerophilic 
vegetation, herbaceous communities in coastal dunes, and 
stands of introduced palm species and banana trees. Ten 
common house geckos, Hemidactylus frenatus DUMERIL 
& BIBRON, 1836, with snout-vent lengths (SVL) = 53.1 + 
1.2 mm, weight = 4.2 + 0.2 g, and seven black-nosed liz- 
ards, Sceloporus melanorhinus BocourT, 1876, SVL = 73 + 
3.6 mm, weight = 18.6 + 2.7 g, were captured on coconut 
palms (Cocos nucifera) and oil palms (Elaeis guineensis); 
five tropical tree lizards, Urosaurus bicarinatus DUMERIL, 
1856, SVL = 44.2 + 1.3 mm, weight = 3.5 + 0.1 g, were cap- 
tured on guanacaste trees (Enterolobium cyclocarpum) and 
fallen trunks. Fifteen terrestrial many-lined whiptail liz- 
ards, Aspidoscelis lineattissimus Core, 1878, SVL = 90.7 + 
1.5 mm, weight = 22.5 + 1.2 g, were captured in open ar- 
eas of tropical deciduous forest and xerophilic vegetation 
where Bursera spp., Cercidium spp., Prosopis spp., Aca- 
cia spp., Opuntia spp. and Stenocereus spp. were the pre- 
dominant plants. Additionally, nine terrestrial red-headed 
spiny lizards, Sceloporus pyrocephalus Cops, 1864, SVL = 
66.6 + 1.5 mm, weight = 14.8 + 0.8 g, were captured in saxi- 
colous habitats and around rock quarries with little or no 
xerophilic vegetation present (Fig. 1). 


Analyses of lizard femoral gland secretions 


Immediately after capture, femoral secretions were ob- 
tained by gently squeezing the femoral pores with forceps. 
Forceps were washed with alcohol and distilled water be- 
tween extractions. Femoral secretions were collected with 
glass inserts and immediately placed inside sterilized glass 
vials that were closed with Teflon-lined stoppers. New 
glass inserts were used for every sample. Vials were stored 
at -20°C until analysis. We also used the same procedure, 
but without collecting secretions, to obtain blank control 
vials that were treated with the same procedure to compare 
with the femoral secretion samples, and to be able to ex- 
clude contaminants from the handling procedure or from 
the environment, and for further examining impurities in 
the solvent. 

Solvent extraction of the solid samples was performed 
adding 200 ul of dichloromethane (Sigma). Samples were 
placed in a vortex for 5 min, then centrifuged for 5 min for 
particle-precipitated elimination (insoluble portion) and 
then stored at -20°C. The solvent-extract phase was collect- 
ed and transferred to appropriate vials for GC-MS analysis. 

The samples were analysed by GC-FID and GC-MS (GC; 
Agilent 7890A, equipped with MSD 5975C and FID detec- 
tors; Agilent Technologies), fitted with a fused-silica cap- 
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illary column (Zebron ZB-FFAP, 30 m length X 0.25 mm 
ID, 0.25-um film thickness; Phenomenex) with helium as 
the carrier gas. 2 ul of each sample dissolved in dichloro- 
methane were injected in splitless mode. Sample injection 
was performed at splitless temperature of 280°C and an 
FID detector temperature of 300°C. Oven temperature was 
programmed to start at 50°C, maintained in isothermal 
for 5 min, then increased to 280°C at a rate of 10°C/min, 
and then isothermal 280°C for 5 min. To analyse using 
both the MS and FID detectors, at the exit of the capillar- 
ity column, a continuous-flow particle separation (50:50) 
was installed. The mass spectrometer was operated at an 
ionization voltage of 70 eV and with scanning between 
m/z 30-500 at 3.9 scans/s. Quantification was performed 
using the relative values of the peak areas in the chroma- 
tograms obtained by using an FID detector. Compounds 
were identified by comparison with the mass spectral li- 
brary NIST/EPA/NIH (ChemStation, Agilent Technolo- 
gies, rev. D.04.00 2014). The synthesized compounds used 
for GC-MS identification consisted of fatty acids (palmit- 
ic and stearic acids) and esters (n-butyl acetate) acquired 
from Sigma-Aldrich. Other synthesized compounds such 
as the esters isobutyl acetate, methyl laurate, isopropyl 
myristate, methyl palmitate, methyl palmitoleate, and me- 
thyl oleate were obtained by chemical synthesis in our lab- 
oratory using the respective chemical precursors. Precur- 
sor compounds were acetic acid, methanol, ethanol, iso- 
propanol, isobutanol, lauric acid, myristic acid, palmitic 
acid, palmitoleic acid, oleic acid, and stearic acid (all from 
Sigma-Aldrich, Merck, or T.J. Baker). Briefly, esters synthe- 
sis was carried out by esterification of their respective al- 
cohols with fatty acids treated with H,SO, and KMnO , at 
130°C for 30 min (CAmpos-Garcia et al. 2018). All the es- 
ter compounds were extracted with dichloromethane sol- 
vent, evaporated, weighed, and dissolved in dichlorometh- 
ane for GC-MS analysis. Additionally, FAME Mix C4-C24 
(catalogue no. 18919, Supelco) were used as standards. 


Statistical analyses 


The proportion of each peak identified in the GC chroma- 
tograms was identified by calculating the relative propor- 
tion of each compound (area under the peak) with respect 
to the total peak area for all chemicals (RAvA-Garcia et al. 
2020). Then, relative areas of the peaks were transformed 
following Aitchison’s formula: [Z = In(Y,/ g(Y,))], where Zi 
is the standardized peak area i for individual j, Y, is the 
peak area i for individual j, and g(Y,) is the geometric mean 
of all peaks for individual j (AITCHISON 1986, DIETEMANN 
et al. 2003). Then, we calculated Euclidean distances be- 
tween every pair of individual samples to produce a re- 
semblance matrix that formed the basis of the analyses. 
We used a single factor permutational multivariate analy- 
sis of variance test PERMANOVA (ANDERSON 2001, MC- 
ARDLE & ANDERSON 2001) based on the Euclidean resem- 
blance matrix using 999 permutations to analyse whether 
the composition of the femoral secretions varied between 
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the five species. All unidentified and identified compounds 
were included in the analysis. Differences between spe- 
cies were further investigated using canonical analysis of 
principal coordinates (CAP) (ANDERSON & WILLIS 2003). 
The software PRIMER V6.1.13 (CLARKE & GORLEY 2006) 
with the PERMANOVA V1.0.3 add-on package (ANDER- 
SON et al. 2008) was used to investigate differences between 
chemical profiles. 


In addition, we assessed if there were statistical differ- 
ences in the relative amount (area under the peak) of each 
compound class among species by fitting negative bino- 
mial GLM analysis (GARDNER et al. 1995). Pairwise com- 
parisons were each performed with a Tukey HSD test. All 
statistical analyses were performed with the interface Rstu- 
dio in R software version 3.3.2 (R Development Core Team 
2013). 


Figure. 1. Species of lizards in their habitats in the tropical Maruata Bay: (A) Hemidactylus frenatus, (B) Sceloporus melanorhinus, 
(C) Urosaurus bicarinatus, (D) Aspidoscelis lineattissimus, and (E) Sceloporus pyrocephalus. 
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Table 1. Chemical compounds identified in femoral gland secretions of arboreal and terrestrial lizard species: Relative proportions 
(mean + SE of the percentage of the total peak areas from FID detector). nd = not detected. RT = retention time and retention index. 


Compounds without standard confirmation (*). 


Arboreal lizards 


Terrestrial lizards 


RT (min) Compound name Rindex H. frenatus S. melanorhinus U. bicarinatus A. lineattissimus S. pyrocephalus 
Esters 

5.2 Isobutyl acetate 1012 0.55£1.32 7.01+9.78 13.27+4.69 9.61+1.65 nd 

6.4 n-Butyl acetate 1074 nd nd 0.22+3.70 0.15+1.05 0.2445.51 

19.3 Methyl laurate 1804 1.46+1.43 31.39+6.33 1.52+5.00 1.28+0.74 1.32+5.00 

22.3 Isopropyl myristate 2027 0.2142.32 0.20+6.19 0.42+1.40 0.3642.13 0.80417.77 

24.5 Methyl palmitate 2208 nd 0.23+3.59 0.97+2.14 0.56+1.03 nd 

24.7 Methyl palmitoleate 2240 2.19+2.69 2.02+8.32 1.23+8.54 1.12+0.65 nd 

27.0 Methyl oleate 2434 nd nd nd 25.64+1.2 nd 
Carboxylic acids 

14.1 Acetic acid 1449 1.35£1;17 0.64+3.70 1.52+7.67 1.18+0.22 4.96+1.59 

31.4 Palmitic acid 2931 39.95+3.04 25.44+6.44 42.62+8.83 25.98+2.50 92.68+5.50 

33.6 Stearic acid 3134 19.27+2.46 11.35+7.13 0.21+1.10 10.89+0.90 nd 
Others 

33.1 t-Squalene* 2865 35.0243.17 21.72+7.00 38.02+8.56 22.93+2.10 nd 

41.3 Dihydro-lanosterol* 3823 nd nd nd 0.29+2.15 nd 

Results We recorded the presence of five compounds in the 


Chemical profiles of lizard species 


We found eight lipophilic compounds in femoral gland se- 
cretions of male lizards H. frenatus (Table 1), which were 
constituted mainly of fatty acid esters corresponding to 
4.4% of total area from all peaks found, three fatty acids 
(acetic, palmitic, and stearic) corresponding to 60.6%, 
and 35% of squalene. The three most abundant chemicals, 
which together comprised > 90% of the total peak areas, 
were palmitic and stearic acid and squalene (Table 1). 

Femoral gland secretions of male S. melanorhinus con- 
tained a mix of nine compounds that could be identified by 
mass fragmentation profiles (Table 1). Overall, these com- 
pounds consisted mainly of esters at 40.9%. These were 
mainly of two types: esters of acetate, such as isobutyl ace- 
tate and esters from fatty acids, such as methyl laurate. The 
three most abundant chemicals, which together comprised 
> 70% of the total peak areas, were methyl laurate, palmitic 
acid and squalene (Table 1). 

In the gland secretions of U. bicarinatus males, ten com- 
pounds including esters at 17.6%, three carboxylic acids at 
43.4%, and others such as squalene at 38% were found. The 
three most abundant chemicals, which together comprised 
> 90% of the total peak areas, were isobutyl acetate, palmit- 
ic acid and squalene (Table 1). 

Femoral gland secretions of A. lineattissimus males con- 
tained a mix of twelve compounds, including esters at 38%, 
carboxylic acids at 38%, and others such as squalene and 
dihydro-lanosterol at 24%. The five most abundant chem- 
icals, which together comprised >90% of the total peak 
areas, were isobutyl acetate, methyl oleate, palmitic acid, 
stearic acid and squalene (Table 1). 
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femoral gland secretions of S. pyrocephalus males, includ- 
ing three esters (2.4%) and two carboxylic acids (97.6%). 
The most abundant chemical was palmitic acid (Table 1). 


Chemical variation between species 

The PERMANOVA analysis based on the resemblance ma- 
trix comparing individuals of each species exhibited signif- 
icant differences in the proportions of compounds between 
the five evaluated species (pseudo Fa = 7:35 P = 0.001). 
The pairwise comparisons revealed significant differences 
in all cases (permutation tests, 2.06 < t < 3.60, p = 0.001 for 
all). The CAP analysis allocated 100% of the chemical pro- 
files to the correct species when using leave-one-out cross- 
validations and m=6 axes (permutation test, 6? = 0.982, 
p = 0.001, Fig. 2). 

With regard to the five species examined, carboxylic 
acids, fatty acid esters and squalene were the most abun- 
dant chemical types, with palmitic acid being the main 
and most abundant compound. We found significant dif- 
ferences in the relative proportions of carboxylic acids be- 
tween our study species (GLM; x’ = 9.63, df = 4, p < 0.001) 
but not in palmitic acid (x? = 1.81, df = 4, p = 0.14). Pair- 
wise comparisons showed that H. frenatus and S. pyro- 
cephalus have lower proportions of carboxylic acids than 
A. lineattissimus and S. melanorhinus (Tukey’s test, p < 0.01 
for all), but U. bicarinatus did not differ significantly from 
any other species (p = 0.15). Squalene, although absent in 
S. pyrocephalus, was the second most abundant compound 
(29%) in the other species. However, we found no signifi- 
cant differences in the relative proportions of squalene be- 
tween species (X? = 2.10, df = 3, p = 0.11). Regarding es- 
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ters, we found significant differences between species (X? = 
9.79, df = 4, p < 0.001), and comparisons revealed that 
S. melanorhinus, U. bicarinatus and A. lineattissimus have 
higher proportions of esters than H. frenatus and S. pyro- 
cephalus (p < 0.01 for all), other comparisons did not yield 
significant differences. 

We found some differences in the chemical profiles be- 
tween arboreal and terrestrial species (Table 1). Arbore- 
al lizards (H. frenatus, S. melanorhinus, and U. bicarina- 
tus) shared eight compounds, whereas the terrestrial liz- 
ards (A. lineattissimus and S. pyrocephalus) shared five 
compounds. Exclusive compounds were found only in 
A. lineattissimus (methyl oleate and dihydro-lanosterol), 
while the remaining species did not have exclusive com- 
pounds. Both groups of species, arboreal and terrestrial, 
shared four compounds; two esters and two carboxylic ac- 
ids (Table 1). 


Discussion 


Our study reports fatty acids, esters, and squalene as the 
main chemical compounds in the femoral gland secretions 
of five tropical lizard species (ALBERTS et al. 1992, MAR- 
TIN et al. 2013a, b, IBÁÑEZ et al. 2017). This composition is 
similar to those recorded from other lizard’s species. For 
example, hexadecanoic acid (palmitic acid) was the main 
fatty acid and the most abundant compound present in 
both terrestrial and arboreal lizards. Other long-chain car- 
boxylic acids, such as octadecanoic acid (stearic acid), were 
found in all studied species except S. pyrocephalus; but pal- 
mitic acid was found in higher proportions than in oth- 
er lizards. Fatty acids of long-chain and very-long-chain 
types have been reported from other lizards from arid and 
desert ecosystems (LOPEZ & MARTIN 2005a, MARTÍN et al. 
2012) and have a functional role as territorial scent markers 
(KHANNOON 2016). Is probable that one or more carboxy- 
lic acids in the secretions of our study lizards facilitate the 
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Figure 2. Canonical analysis of principal coordinates (CAP), 
showing classification of chemicals identified in femoral secre- 
tions of three arboreal H. frenatus, S. melanorhinus, and U. bi- 
carinatus; and two terrestrial A. lineattissimus and S. pyrocephalus 
lizard species. 


development of similar social contexts, but this is as yet an 
information gap that needs to be addressed. 

Some lizards inhabiting high mountain areas with low 
temperatures and high humidity frequently have greater 
relative proportions of esters in their femoral gland secre- 
tions than lizards inhabiting low-altitude areas (GABIROT 
et al. 2012, MARTÍN et al. 2017). This ester variation is ex- 
plained as an adaptation to maximize efficacy of substrate 
scent marks under different microclimatic conditions 
(MARTÍN et al. 2015, BAECKENS et al. 2018). In our study, 
lizards from tropical climates with high temperatures and 
humidity produced relatively higher amounts of esters (5- 
42%) than the high-mountain lizard Psammodromus algi- 
rus (5.5% of esters) from Mediterranean climates (MARTÍN 
et al. 2017). Possibly based on its chemical properties, esters 
are more volatile and prone to oxidation than carboxylic 
acids in environments with tropical climates. Therefore, it 
is possible, but has not as yet been conclusively demon- 
strated, that more humid habitats are conducive to higher 
proportions of esters due to accelerated hydrolysis. 

Although in relatively low proportions, squalene is a 
common compound in femoral secretions of many lizard 
species (WELDON et al. 2008). This chemical is considered 
an antioxidant that will enhance the signalling function of 
other compounds (MARTÍN & LÓPEZ 2014). As squalene 
might stabilize other lipid compounds by limiting oxida- 
tion (GABIROT et al. 2008), its abundant presence in four of 
the five studied species may be an adaptation to the usual- 
ly high humidity of their tropical environment. Squalene’s 
role as a pheromone has only been tested in fossorial am- 
phisbaenians, and it has been considered an important 
component in antagonistic and sex-discrimination behav- 
iours (LÓPEZ & MARTÍN 2005b). 

Although diverse steroids are frequently present in the 
femoral secretions of lizards (WELDON et al. 2008), in our 
study only one steroid, the lanost-8-en-3ß-ol (dihydro- 
lanosterol), was found in one terrestrial species (A. lineat- 
tissimus). Cholesterol and some of its derivatives are the 
most important steroids recorded in almost all lizard spe- 
cies (MARTÍN et al. 2015). However, there are other ster- 
oids of plant origin (phytosterols), one example of these 
is lanosterol, abundant in the femoral gland secretions of 
the green iguana, Iguana iguana, suggesting a possible re- 
lationship between diet and the chemical characteristics 
of gland secretions (ALBERTS et al. 1992, BAECKENS et al. 
2017). Similar to lanosterol, dihydro-lanosterol, present in 
A. lineattissimus, is a phytosterol that could possibly have 
some relationship with the food spectrum of this species. 

In conclusion, we identified a significant diversity of li- 
pophilic compounds in the femoral gland secretions of ar- 
boreal and terrestrial species of tropical lizards and clear 
interspecific differences in the fatty acids and esters these 
contain. Some of the compounds reported in this study 
have been described as semiochemicals with signalling 
function (pheromones) involved in the sexual communi- 
cation of some lizard species, but their behavioural role in 
tropical lizards remains to be confirmed through future re- 
search. 
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